Probing the Nuclear Equation of State in the quasi-elastic 
nucleus-nucleus scattering 
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Large-angle elastic scattering of a-particle and strongly-bound light nuclei at a few tens MeV/nucleon 
has shown the pattern of rainbow scattering. This interesting process was shown to involve a significant 
overlap of the two colliding nuclei, with the total nuclear density well above the saturation density of 
normal nuclear matter (NM) . For a microscopic calculation of the nucleus-nucleus potential within the 
folding model, we have developed a density dependent nucleon-nucleon (NN) interaction based on the 
G-matrix interaction M3Y. Our folding analysis of the refractive ^He, ^^C, and ^^O elastic scattering 
shows consistently that the NM incompressibility K should be around 250 MeV which implies a rather 
soft nuclear Equation of State (EOS). To probe the symmetry part of the nuclear EOS, we have used 
the isovector coupling to link the isospin dependence of the proton optical potential to the cross section 
of {p, n) charge-exchange reactions exciting the isobaric analog states in nuclei of different mass regions. 
With the isospin dependence of the NN interaction fine tuned to reproduce the charge exchange data, 
a realistic estimate of the NM symmetry energy has been made. 



1 What is the nuclear rainbow? 

The atmospheric rainbow is observed in 
nature whenever there are water droplets illu- 
minated by the sun light. It can be seen dur- 
ing the rain with the sunshine not completely 
covered by the clouds or from a fountain when 
the sunlight enters from behind the point of 
observation. Beside the fascinating effect of 
color splitting in the rainbow caused by the 
dependence of the refraction index on the light 
wavelength, a physically more interesting ef- 
fect is the increased light intensity around the 
rainbow angle Or and the shadow region lying 
beyond G/j which are results of a particular 
refraction - reflection sequence. 

Although Descartes has successfully ex- 
plained the origin of the atmospheric rain- 
bow based on simple geometrical ray optics 
in 1637, neither he nor Newton and Young 
(several decades later) could explain the fine 
structure of the rainbow seen as the supernu- 
meraries (the faint bows located just below 
the primary bow). The first complete math- 
ematical description of the atmospheric rain- 
bow was given in 1838 by Airy, and the oscil- 
lation of the light intensity near the rainbow 
angle is now known as Airy oscillation which 



gives rise to the supernumeraries. As nuclei 
are approximately spherical objects (like wa- 
ter drops) having wave properties, they can 
be refracted or undergo interference in the 
nucleus-nucleus scattering, like the refraction 
of sunlight by the water drops. As a result, 
one may observe phenomena like rainbow in 
the nuclear scattering if the conditions are 
right. Indeed, the rainbow pattern was clearly 
observed in large-angle elastic scattering of 
a-particles and some strongly-bound light 
nuclei at few tens MeV/nucleon, where sig- 
natures of the Airy oscillation pattern has 
been identified. An important feature of the 
nuclear rainbow is that it allows us to probe 
the nucleus-nucleus interaction at small dis- 
tances [1] thanks to a weak absorption in the 
nucleus-nucleus system. 

2 From nuclear rainbow to the equa- 
tion of state for cold nuclear matter 

For a microscopic description of the elas- 
tic nucleus-nucleus scattering, the folding 
model analysis is usually performed, where 
the (real) nucleus-nucleus optical potential is 
calculated Hartree-Fock (HF) potential 
of the dinuclear system [2,3] using an effective 



nucleon-nucleon (NN) interaction [4] 

^ = J2i< ^J>d|u > + < ij\vEx\ji >], (1) 

where \i > and \j > are the single-particle 
wave functions of nucleons in the two collid- 
ing nuclei Ai and A2, respectively; vn and v ex 
are the direct and exchange parts of the ef- 
fective NN interaction. It turns out that the 
density dependence of the effective NN inter- 
action can be accurately tested in the fold- 
ing analysis of the refractive elastic a-nucleus 
or nucleus-nucleus scattering data [2,5]. For 
this purpose, a phenomenological density de- 
pendence was first introduced to the M3Y in- 
teraction based on the G-matrix elements of 
the Reid and Paris NN potentials [6], to re- 
produce to the saturation properties of cold 
nuclear matter (NM) in the HF scheme [4]. 
From the HF results for the NM energy (per 
nucleon) plotted in Fig. 1, one can see that 
different sets of the density dependence give 
values of the NM incompressibility K rang- 
ing from 170 to above 500 MeV. Since the K 
value is a key input in the NM equation of 
state (EOS), a test of the density dependence 
of the NN interaction is also an indirect test 
of the nuclear EOS. In the elastic channel, the 
two colliding nuclei remain in the ground state 
(g.s.) even when they overlap strongly at small 
impact parameters, since any density defor- 
mation or rearrangement directly projects the 
system out of the elastic channel. Therefore, 
the total density p of the two overlapping nu- 
clei which enters Eq. (1) must be taken as the 
sum of the two g.s. densities and the total den- 
sity for a projectile overlapping a target nu- 
cleus may reach as much as twice the normal 
NM density po [1;2]. Note that the shorter the 
impact parameter the higher the overlap den- 
sity and the more the nucleus-nucleus poten- 
tial becomes sensitive to the density depen- 
dence of the NN interaction. At small impact 
parameters (or low partial waves) the survival 
probability of the elastic wave is usually less 
than 1% and requires, therefore, a very precise 
measurement of elastic events at large scatter- 
ing angles. These large-angle data points give 
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Fig. 1. EOS for cold NM given by the HF calcu- 
lation using different density dependent NN in- 
teractions (which give different values of the NM 
incompressibility K). 

us the most vital information about the reli- 
ability of the effective NN interaction used in 
the folding calculation [8] . From results of the 
folding analysis of the elastic ^^O+^^O scat- 
tering data at 350 MeV [7] plotted in Fig. 2 
one can clearly see which is the most appropri- 
ate density dependence of the NN interaction 
(in the corresponding K value). Together with 
the results obtained earlier for elastic scatter- 
ing of a particles and other light projectiles 
[8], we conclude that the most realistic value 
for the NM incompressibility K is around 230 
- 270 MeV which corresponds to a rather soft 
EOS. We stress that the broad bump seen in 
Fig. 2 at 6c.m. ~ 50° has been specified as the 
first Airy maximum [9] , and the observed rain- 
bow pattern is quite sensitive to the ^^O+^^O 
optical potential at small distances. It can be 
shown for the considered ^^0-1-^^0 system (see 
Fig. 3) that the large-angle data are sensitive 
to the impact parameters as small as 2 fm. 
At such small internuclear distances, the total 
overlap density of the system reaches up to 2po 
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Fig. 2. Folding model description of the elastic 
160+160 scattering data at E'lab = 350 MeV [7 
using the same density dependent NN interac- 
tions as those in Fig. 1. The best-fit interaction 
is CDM3Y6 which gives K ^ 252 MeV. 

[8]. Therefore, the elastic refractive nucleus- 
nucleus scattering data like those measured 
for the ^^0-|-^^0 system at 350 MeV provide a 
very good data base for probing the high over- 
lap density in the elastic channel. 

3 Probing EOS of the asymmetric NM 
via charge-exchange reaction 

The nuclear EOS presented in Fig. 1 was 
obtained in a HP calculation for symmetric 
NM, i.e., with equal proton and neutron den- 
sities. In reality, the NM that exists in the 
neutron stars is highly asymmetric, with the 
neutrons outnumbering protons by factor of 
about 2 in the crust of a neutron star. There- 
fore, the knowledge about the symmetry part 
of the EOS is vital for the understanding of 
the dynamics of supernovae explosion and the 
formation of neutron stars [10,11]. The sym- 
metry part of the nuclear EOS is determined 
essentially by the NM symmetry energy S{p) 
(defined in terms of a Taylor series expansion 
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Fig. 3. The same description as in Fig. 2 given 
by the CDM3Y6 interaction, but with different 
cutoff values of the lowest partial wave Lmin (and 
the corresponding impact parameters Rmin)- 

of the NM binding energy B = E/A) as 

S(p, 5) = Bip, 0) + Sip)6' + 0(5') + ... (2) 

where 6 = {pn — Pp) /p is the neutron-proton 
asymmetry parameter. The contribution of 
0{S^) and higher-order terms in Eq. (2), i.e., 
the deviation from the parabolic law was 
proven to be negligible [12]. The NM symme- 
try energy determined at the NM saturation 
density, -Egym = S{po) with po ~ 0.17 fm~^, is 
widely known in the literature as the symme- 
try energy or symmetry coefficient. Although 
numerous nuclear many-body calculations 
have predicted -Esym to be around 30 MeV, 
a direct experimental determination of -Egym 
still remains a challenging task. One needs, 
therefore, to relate -Egym to some experimen- 
tally inferrable quantity like the neutron skin 
in neutron- rich nuclei [13] or the fragmenta- 
tion data of heavy-ion (HI) collisions involv- 
ing N Z nuclei [14]. Within the frame of 
any microscopic model for asymmetric NM, 
the symmetry energy depends strongly on 
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the isospin dependence of the NN interaction 
used therein [12]. Therefore, the Esym value 
can be indirectly tested in a charge exchange 
(isospin-fiip) reaction which has been known 
for decades as a good probe of the isospin de- 
pendence of the effective NN interaction [15]. 
Although the isospin dependence of the nu- 
clear optical potential (OP), known by now 
as Lane potential [16], has been studied since 
a long time, there has been a considerable in- 
terest recently in studying the isospin depen- 
dence of the OP in the quasi-elastic scattering 
reactions measured with unstable neutron- 
rich beams. Based on the isospin symmetry, 
the nucleus-nucleus OP can be written in 
terms of an isovector coupling [16] as 

U{R) = Uo{R)+AUi{R)^-^, (3) 

where t is the isospin of the projectile a and 
T is that of the target A. While the rela- 
tive contribution by the Lane potential Ui 
to the elastic {p,p) cross section is small and 
amounts only to a few percent for a neutron- 
rich target, it determines entirely the (Fermi- 
type) AJ'^ = 0"^ transition strength of the 
(p, n) reaction leading to an isobaric analog 
state (IAS). Therefore, the {p,n) reaction has 
been so far the main tool in studying the 
isospin dependence of the proton-nucleus OP. 
Since this isospin dependence should be bet- 
ter tested in the charge exchange reactions 
induced by the neutron-rich beams, we con- 
sider in the present work the p(^He,^Li*)n 
reaction measured by Cortina- Gil et al. [17] 
with the secondary ^He beam at E'lab = 41.6^4 
MeV. Given a large neutron-proton asymme- 
try (5 = 1/3) of the unstable ^He nucleus, 
the measured p(^He,^Li*)n cross section for 
the transition connecting the ground state 
of ^He (T = = 1) and its isobaric ana- 
log partner (T = l,Tj = 0, J'" = 0"*" excited 
state of ^Li at 3.563 MeV) is indeed a good 
probe of the isovector coupling in the ^He+p 
system. To link the Lane potential Ui to the 
isospin dependence of the NN interaction, 
we have used the folding model [18] to cal- 
culate Uq and Ui using the explicit proton 
and neutron g.s. densities of ®He and the 



CDM3Y6 density- and isospin dependent NN 
interaction [8]. The only nuclear structure in- 
put is the ^Hcg.s. density and we have used 
the microscopic density given by the cluster- 
orbital shell model approximation (COSMA) 
[19]. The p(^He,^Li*)n cross sections given 
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Fig. 4. CC results for the charge exchange 
p(^He,^Li*)n cross section at E'lab = 41.6^ 
MeV in comparison with the data measured by 
Cortina-Gil et al. [17]. 

by the coupled-channel (CC) calculation us- 
ing a charge-exchange form factor based on 
the Lane potential Ui were found to have 
a shape very close to that of the measured 
angular distribution (see Fig. 4). Since the 
complex strength of the form factor was fixed 
by the folding model analysis of the elastic 
^He+p scattering [20], the CC description 
of the p(^He,^Li*)n data could be improved 
only by fine tuning the strength Ci of the 
isovector part of the density dependence of 
the CDM3Y6 interaction [12]. One can see 
that the best fit is achieved when Ci is about 
10% stronger than the isoscalar strength Cq. 
To make a more definitive conclusion on the 
EOS of asymmetric NM, we found it neces- 
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Fig. 5. CC results for the charge exchange 
^^Ca(p, n)^^Sc cross section at the incident pro- 
ton energy of 45 MeV in comparison with the 
data measured by Doering et al. [15]. 

sary to make a systematic folding analysis of 
the charge exchange {p, n) reactions measured 
with targets in different mass regions. There- 
fore, the same density- and isospin dependent 
NN interaction has been used to construct the 
charge exchange form factors for the (p, n) 
reactions measured at the incident proton 
energies of 35 and 45 MeV with the targets 
^«Ca, 90Zr, i20Sn and ^ospb [15]. Although the 
neutron-proton asymmetry 5 of these nuclei 
is smaller than that of unstable ^He nucleus, 
the complex proton OP for these nuclei has 
been studied over the years [21] and, hence, it 
allows us to unambiguously probe the isovec- 
tor part of the OP via the (p, n) reaction. By 
using a most appropriate complex proton OP 
which produces not only the elastic scattering 
data but also the polarization data and the 
experimental total reaction cross section for 
each of the considered targets, we have come 
up with about the same accurate description 
of the (p, n) reactions leading to the excitation 
of IAS in ^^Ca, ^^Zr, ^^ogn and ^o^Pb as that 



presented above for the unstable ^He. Our CC 
results for the charge exchange ^^Ca(p, n)^^Sc 
cross section at 45 MeV and the data mea- 
sured by Doering et al. [15] are shown in 
Fig. 5 as an illustration of the success of our 
approach. The best-fit isovector strength Ci 
of the CDM3Y6 interaction was adjusted in 
each case not only to reproduce the measured 
angular distribution for the considered (p, n) 
reaction but also to obtain the measured total 
{p, n) cross section in our CC calculation. For 
example, the Ci values used to calculate the 
^^Ca(p, n)^^Sc cross section shown in Fig. 5 
was constrained by a CC calculation giving 
the experimental total ^^Ca(p, n)'^^Sc cross 
section of 8.4 ± 1.0 mb [15]. If one takes into 
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Fig. 6. Density dependence of the NM symmetry 
energy S{p) predicted by the HF formalism [12] 
using the same isovector strengths Ci as those 
used in Fig. 4 and the (empirical) neutron-skin 
[13] and HI fragmentation [14] data. 

account only the uncertainty of the measured 
(p, n) angular distribution, the range of ac- 
ceptance for Ci values becomes significantly 
larger. Given the results of our CC calcu- 
lation of the (p, n) reactions for '^^Ca, ^°Zr, 



^20Sn and ^nspb [22], we conclude that the 
best-fit isovector strength Ci of the CDM3Y6 
interaction is slightly larger than the isoscalar 
strength Cq, with about the same uncertainty 
as that found in the analysis of p(^He,^Li*)n 
reaction. 

With the isovector strength of the 
CDM3Y6 interaction now well tested, we have 
further performed the HF calculation [12] of 
asymmetric NM using this same isospin- and 
density dependent interaction. The density 
dependence of the NM symmetry energy S{p) 
obtained with the same isovector strengths 
Ci as those used in Fig. 4 is shown in Fig. 6, 
and one can deduce easily £^sym ~ 31 ± 1 
MeV from our HF results. This result is quite 
complementary to the structure studies which 
relate the Egym value to the neutron skin, a 
method first suggested by Brown [23] . If one 
adopts a neutron-skin AR Ri 0.1 — 0.2 fm for 
^°^Pb then a systematics based on the mean- 
field calculations [13] gives -Egym ~ 27 — 31 
McV (which is plotted in Fig. 6). Our result 
is also complementary to the recent studies 
of HI fragmentation based on the antisym- 
metrized molecular dynamics [14] which ob- 
tained S{p ^ 0.08 fm-3) ^18-22 MeV at 
a finite temperature around 3 MeV. If we ne- 
glect the temperature dependence of S{p) at 
low NM densities, this value turns out to agree 
well with our HF result for the low-density 
part of S{p) as shown in Fig. 6. 

In conclusion, the charge exchange reac- 
tion like (p, n) or (^He,t) should be a very 
powerful tool in studying not only the nuclear 
structure but also the isospin aspects of the 
EOS for asymmetric NM. 
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